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A single-shot Toffoli, or controlled-controlled-NOT, gate is desirable for classical and quantum 
information processing. The Toffoli gate alone is universal for reversible computing and, accompa¬ 
nied by the Hadamard gate, forms a universal gate set for quantum computing. The Toffoli gate 
is also a key ingredient for (non-topological) quantum error correction. Currently Toffoli gates are 
achieved by decomposing into sequentially implemented single- and two-qubit gates, which requires 
much longer times and yields lower overall fidelities compared to a single-shot implementation. 

We develop a quantum-control procedure to construct a single-shot Toffoli gate for three nearest- 
neighbor-coupled superconducting transmon systems such that the fidelity is 99.9% and is as fast 
as an entangling two-qubit gate under the same realistic conditions. The gate is achieved by a 
non-greedy quantum control procedure using our enhanced version of the Differential Evolution 
algorithm. 

PACS numbers: 03.67.Lx, 85.25.Cp, 42.50.Ex 


Scalable quantum computing [1, 2] requires a set of 
high-fidelity universal quantum gates with which to con¬ 
struct the circuit [2-4]. Experimental progress towards 
a high-fidelity universal set of gates comprising single- 
and two-qubit operations has been impressive, exceeding 
99.9% for single-qubit gates and 99% for an entangling 
two-qubit gate [3], but an outstanding problem is that 
(non-topological) quantum error correcting codes require 
a gate acting on at least three qubits [5, 6], with the Tof- 
foli gate [7-9] being optimal. The Toffoli gate is also a 
key component for reversible arithmetic operations, such 
as the modular exponentiation, which is a necessary step 
in Shor’s factoring algorithm [10]. 

The quantum Toffoli gate is to effect a three-qubit 
controlled-controlled-NOT (CCNOT) gate, which means 
that the third qubit is flipped only if the first two qubits 
are in the |1) state and not flipped otherwise. Thus far 
Toffoli gates are achieved by decomposing into single- and 
two-qubit gates with resultant fidelities limited to 81% 
in a post-selected photonic circuit [11], 71% in an ion- 
trap system [7], 68.5% in a three-qubit circuit QED sys¬ 
tem [8], and 78% in a four-qubit circuit QED system [6]. 
We here introduce a non-greedy quantum-control ap¬ 
proach for directly constructing Toffoli gates based on 
an enhanced version of the Differential Evolution (DE) 
algorithm [12, 13]. We show that our scheme applied 
to the three nearest-neighbor-coupled superconducting 
transmon systems should produce a Toffoli gate oper¬ 
ating with 99.9% fidelity and operating as fast as an en¬ 
tangling two-qubit gate under the same conditions. As 
our quantum-control-based approach [14] to realizing the 
Toffoli gate does not resort to decomposition, a fast Tof¬ 


foli gate enables error-correction with high fidelity under 
this scheme. An additional valuable benefit of realizing 
CCNOT directly is that the Hadamard (H) and CCNOT 
together make a universal gate set [15] with signihcant 
advantages over the oft-studied H, tt/S gate and CNOT 
universal set [2]. 

Superconducting circuits offer a promising medium for 
realizing a high-fidelity CCNOT gate based on quan¬ 
tum control of three nearest-neighbor-coupled supercon¬ 
ducting artificial atoms [3]. Our approach is to vary 
the energy levels for each of three individual supercon¬ 
ducting atoms using time-dependent control electron¬ 
ics, which conveniently do not require additional mi¬ 
crowave control [16]. A similar strategy has recently been 
employed successfully to design two-qubit controlled-Z 
gates, for which optimal pulses are found via a greedy 
algorithm [17]. We, however, have observed that ex¬ 
isting optimization algorithms (including greedy algo¬ 
rithms) are insufficient to generate an optimal pulse 
for high-fidelity Toffoli gates, and, therefore, developed 
a non-greedy optimization scheme, referred to here as 
Subspace-Selective Self-Adaptive DE or SuSSADE. 

We consider a linear chain of three nearest-neighbor- 
coupled superconducting artificial atoms, realized as 
transmons [3] with distinct locations labeled k = 1,2,3. 
The transmons have non-degenerate discrete energy lev¬ 
els, labeled {|j)fc}, with j = 0 for the ground state. The 
energies are anharmonically spaced, with this spacing al¬ 
lowed to be dependent on the specihc transmon. Whereas 
superconducting atoms contain many energy levels, we 
truncate all energy levels for j > 3 for each transmon as 
a CCNOT operates on at most three excitations. 
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Although the TofToli gate acts on three qubits as per 
definition, our quantum-control procedure operates on 
the first four levels of each transmon. The Hamiltonian 
that generates Toffoli acts on the 4^-dimensional Hilbert 
space with energy basis {\j G {0,1,2,3})®3}. We 

follow the standard practice of specifying transmon en¬ 
ergy levels instead as frequencies with these atomic fre¬ 
quencies shifted by frequency of a rotating-frame basis 
transformation: the shifted frequency of the trans¬ 
mon is Afc, and the anharmonicity of the j**' level of 
the transmon is rjjk. Therefore, the energy of the fc**' 
transmon’s level at time t is h{jAk{t) — rjjk). 

Nearest-neighbor transmons couple via an XY interac¬ 
tion with coupling strength between the fc*'' and (fc-|-1)**' 
transmons denoted by g^- The three-transmon Hamilto¬ 
nian is thus [16] 

- Vjk) \j)k O'lfe 

k—1j—0 
2 

+ ^ (AfeAfc_|_i-I-yfeYfc_|_i), (1) 

k=l ^ 

for coupling operators 

3 

Xk='^'/j\j - l)fc (j|fe +hc, 

1=1 

3 

Yk=-Yl \/^\j - l)fc {j\k + he, (2) 

1=1 

which are higher-dimensional generalizations of Pauli 
spin matrices [16, 18], and he denotes Hermitian con¬ 
jugate. 

Here we employ Hamiltonian evolution to realize the 
CCZ gate, which effects ajO) +/311) >->■ ajO) — /S]!) on the 
third qubit only if the first two qubits are jl). The CC- 
NOT and CCZ operations are equivalent under the local 
transformation CCNOT= [ll(g)l 0 H] CCZ [l(g)l 0 H] 
(similar to the equivalence between two-qubit CNOT and 
CZ gates), with H straightforward to implement as a 
fast single-qubit operation [19, 20]. The CCZ gate is 
achieved by varying of each superconducting atom 
over duration 0 with resultant Hamiltonian-generated 
time-ordered (T) evolution operator 

17(0) = Texp i J 77(r)dr^ . (3) 

Whereas our approach enables generating any desir¬ 
able pulse shape for A^, here we consider two types of 
pulses: piecewise-constant and piecewise-error-function. 
These time-dependent control pulses are constrained 
within the frequency-range of a superconducting trans¬ 
mon system. We employ the less computationally- 
expensive piecewise-constant function to demonstrate the 


existence of an optimal pulse for high-fidelity Toffoli gate. 
However, the control electronics for superconducting sys¬ 
tems is only capable of generating smooth pulses, which 
motivated us to consider a realistic case for which the 
control parameters are connected together via smooth 
error functions [16, 21]. We show that the gate fidelity 
does not depend on the pulse shape, and only depends 
on the number of control parameters. Therefore, without 
any loss of generality, we choose the less computationally 
expensive piecewise-constant control function to analyze 
the fidelity of the designed gate against other parameters. 

The Hamiltonian evolution (3) describes the unitary 
dynamics of the system in the absence of decoherence. 
Decoherence is incorporated by treating each atom as 
a damped harmonic oscillator characterised by ampli¬ 
tude and scattering induced phase-damping rates for 
each oscillator. The corresponding timescales are relax¬ 
ation time Ti and dephasing time T 2 , analogous to the 
rates employed for two-level systems [2, 22]. We here 
assume T := Ti = T 2 , which is valid for frequency- 
tunable transmons [16]. These decohering processes 
modify the unitary evolution (3) to a completely-positive 
trace-preserving map i£’(0), which is decomposable into 
an operator sum as discussed in the Supplementary Ma¬ 
terial [23]. 

A high-fidelity quantum gate is usually designed by 
determining an optimal control pulse for each fre¬ 
quency Afc neglecting open-system effects such as de¬ 
coherence. Performance is assessed for the unitary 
evolution (3) projected to the computational subspace: 
17^(0) := The standard figure of merit 

for performance of U is the “intrinsic fidelity” (fi¬ 
delity neglecting decoherence) with respect to the ideal 
gate, in this case CCZ, so the intrinsic fidelity is [13] 
F = I Tr^CCZl' 17^^(0)) with = 1 if 17^(0) = 


CCZ and 0 < < 1 otherwise. After determining con¬ 

trol pulses that maximize F, decoherence is then incor¬ 
porated into the calculation to assess the performance 
under open-system conditions [16]. 

In the presence of decoherence, the efficacy of the non¬ 
unitary evolution compared to the target gate is quan¬ 
tified by the average state fidelity F, which is calcu¬ 
lated as follows. For {['^fc) S the set of three- 

transmon computational basis states, the non-unitary ex¬ 
tension of the unitary evolution (3) transforms a pure 
computational basis state IV’fe) (V’fel into a mixed state 
As each basis state \'tpk) {4’k\ remains invari¬ 
ant under an ideal CCZ gate, average state fidelity 

= I Hk '\J\{i’k\ \i’k)\ quantifies the efficacy of 

a quantum gate in the presence of intrinsic as well as 
decoherence-induced noise for a given optimal pulse. 
Whereas F « 99.9% is considered to be a threshold for 
topological (surface-code) fault-tolerance for single- and 
two-qubit gates [3], our approach achieves this fidelity 
even for the three-qubit CCZ gate subject to realistic 
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constraints of the control pulses. In this work, unless 
otherwise stated, the average state fidelity is referred to 
as fidelity. 

The strategy for controlling the evolution (3) is to 
vary the frequencies so that energy levels approach each 
other but then avoid degeneracies, known as avoided level 
crossings. These avoided crossings mix energy-level pop¬ 
ulations and dynamical phases together. This avoided- 
crossing effect enables shaping the evolution toward the 
final time-evolution operator objective, which is obtained 
by maximizing T. 

Optimal pulse shapes for each are obtained by dis¬ 
cretizing the time duration 0 into iV constant intervals of 
duration At := 0/A, and the control-problem parame¬ 
ter space is spanned by the set of variables { Afe(^AT); £ = 
1,..., A} for each k. These control-points are then con¬ 
nected via step functions or error functions to construct 
the pulse shapes as we described earlier. The CCZ opti¬ 
mization problem is non-convex with a 3A-dimensional 
parameter space corresponding to A parameters for each 
of the three frequencies Ai^ 2 , 3 - For a fixed At, therefore, 
the dimension of the parameter space increases linearly 
with the total time duration 0, which influences which 
optimization methods work and which do not. 

We devise a quantum-control procedure that designs 
an optimal pulse for a Toffoli gate, which operates as fast 
as a two-qubit gate [16] with a target intrinsic fidelity 
of 0.9999. We first use the existing optimization algo¬ 
rithms, namely quasi-Newton approach, which employs 
the Broyden-Fletcher-Goldfarb-Shanno (BFGS) approxi¬ 
mation of the Hessian [24-28] , simplex methods [29] , sev¬ 
eral versions of particle swarm optimization [30, 31], and 
differential evolution (DE) algorithms [13, 32]; however, 
they all failed to reach the target intrinsic fidelity under 
the time-constraint of the fast Toffoli gate. Therefore, we 
construct a new optimization algorithm here to realize a 
Toffoli gate that reaches our target. 

Of these optimization approaches, DE yielded the best 
fidelity but failed to reach the target due to the well 
known problems of searching high-dimensional parameter 
spaces [33]. This drawback motivated us to enhance DE 
for such high-dimensional problems by instead breeding 
over randomly selected low-dimensional subspaces, hence 
our name Subspace-Selective Self-Adaptive DE (SuS- 
SADE) algorithm (see Supplementary Material [23]). 
One of our objectives is to demonstrate (numerically) 
the capability of SuSSADE with respect to finding a so¬ 
lution equally successfully regardless of parameter-space 
dimension within the regime that is relevant for current 
superconducting experiments. 

To understand our enhancement, we first briefly review 
standard DE [32]. DE cooperatively evolves a collection 
of trial solutions, called chromosomes, towards an opti¬ 
mal solution. Ghromosomes are labeled by their location 
in the parameter space, and optimization is thus a search 
for the best chromosome in this space. Evolution from 


one generation (i.e., chromosomes for one iteration step) 
to the next is achieved by breeding each chromosome 
with three other randomly chosen chromosomes from the 
same generation. Breeding yields a single daughter chro¬ 
mosome, and only the fittest of the original and daughter 
chromosome survives. This breeding-and-survival proce¬ 
dure continues until either a chromosome reaches the req¬ 
uisite or the number of generations reaches a specified 
upper bound. 

Whereas standard DE breeds chromosomes randomly 
selected from the entire space, our SuSSADE algorithm 
is much faster due to breeding being restricted some 
of the time to a subset of chromosomes drawn from a 
low-dimensional subspace, i.e., some fixed parameters 
and some variable parameters. Our algorithm randomly 
switches breeding between the subspace and the whole 
space according to the value of an input switch parame¬ 
ter S S [0,1] such that a uniformly distributed random 
number Vj € [0,1] at generation j restricts breeding to 
the subspace if Vj < S and breeds in the whole space 
otherwise. 

In the extreme case of restricting to one-dimensional 
subspaces, chromosomes can breed only if all but one 
of the parameters are the same. We refer to this one¬ 
dimensional extreme case as IDSuSSADE. Henceforth we 
use only IDSuSSADE as it works well with S = 0.14 for 
designing the Toffoli gate. 

Here we present two types of pulses that achieve the 
target intrinsic fidelity, and we explore how the perfor¬ 
mance of piecewise constant pulses vary with respect to 
the total gate time, coupling strength and decoherence- 
induced noise. The success of our quantum control pro¬ 
cedure corresponds to a target intrinsic fidelity of 0.9999 
and a timescale comparable to a two-qubit gate [16]. 

Figure 1 shows both piecewise-constant as well as 
piecewise-error-function pulses that achieve the target 
intrinsic fidelity obtained by optimizing all the param¬ 
eters within the experimental constraints. The GGZ 
gate corresponding to Figure 1 requires a total gate time 
of 26 ns given a coupling strength of g = 30 MHz [16, 17]. 
Gomparing the intrinsic fidelities of Figure 1(a) and Fig¬ 
ure 1(b) shows that the target fidelity does not depend 
on the shape of the pulse; rather it depends on the num¬ 
ber of control parameters. In what follows, therefore, 
we consider only the piecewise-constant pulse shapes as 
these are computationally less expensive to handle and 
also do not compromise the generality of our results. We 
show in Figure 2 how the (maximized) intrinsic fidelity 
changes when the parameters g and 0 are varied within 
a range commensurate with currently available super¬ 
conducting circuits [3]. Figure 2(a) gives the intrinsic 
fidelity as a function of total gate time for various cou¬ 
pling strengths g. Figure 2(b) shows that the total gate 
time changes linearly with the coupling g, in order to 
achieve a given fidelity. Finally, we consider the effect 
of decoherence on the approximate CCZ gate obtained 


4 


A 


1 




FIG. 1: (color online) Optimal pulse shapes for the Toffoli 
gate given as frequency detunings (for the i = 1,2,3) of 
the superconducting atoms, corresponding to A) piecewise- 
constant and B) error-function-based pulse profiles, as a func¬ 
tion of time r with constant step-size time interval Ar = 1 
ns and with T = 0.9999 and g = 30 MHz. The black dots on 
both plots show the control parameters used to optimize the 
shape of the pulses. 



FIG. 2: (color online) A) Intrinsic fidelity vs. total gate 
time 0 for various coupling strengths g and B) 1/0 vs g for 
J- — 0.999. The O, A, o, and □ denote the actual numerical 
computations using IDSuSSADE, and solid lines depict cubic- 
fit curves. 


by optimal pulses shown in Figure 1(a), and compute 
the fidelity T. Amplitude-damping and phase-damping 
rates and are treated as the dominant forms 

of decoherence. For fast gates with 0 ^ T, but with de¬ 
coherence more significant than intrinsic errors, an order- 
of-magnitude estimate yields 1— ^ ^ 0/T [16], which is 
consistent with the numerically evaluated plot of .F vs T 
in Figure 3. We have employed our quantum control 
procedure to determine the optimal pulse (Figure 1) for 
a high-fidelity single-shot three-qubit Toffoli gate. We 
computed a smooth pulse (Figure 1(b)), for which the 
control parameters are separated by 1 ns and connected 
via error functions, thereby ensuring that the pulse is 
compatible with the power and bandwidth specifications 
of standard control electronics. Applying our approach 
to the three nearest-neighbor-coupled superconducting 
transmon systems produces a fast and high-fidelity Tof¬ 
foli gate in 26 ns, which matches the timescale for the 
two-qubit avoided-crossing-based CZ gate. 



T(gs) 


FIG. 3: (color online) The fidelity F is plotted against the 
coherence times T. We assnme T = Ti = T 2 , with Ti and 
T 2 the relaxation time and dephasing time of each transmon 
respectively. This assumption is valid for tnnable transmons. 
Each O denotes an actual numerical result obtained from the 
decoherence calculation. The solid line depicts the cubic-fit 
curve. 


The longer total gate time, with a fixed Ar, gener¬ 
ates a higher-dimensional parameter space for the opti¬ 
mization algorithm. The monotonically increasing op¬ 
timized intrinsic fidelity in Figure 2(a) with increasing 
total time thus demonstrates the capability of our algo¬ 
rithm for a parameter regime relevant to superconducting 
experiments, for which alternative algorithms fail. The 
linear relationship between 1/0 and g demonstrates that 
faster gate speed requires higher coupling. This relation 
is a characteristic signature for avoided-crossing-based 
gates (as also obtained for avoided-crossing-based two- 
qubit gates [16]) assuming the corresponding optimiza¬ 
tion algorithm is capable of finding the optimal solution 
regardless of the parameter space dimension. 

The effect of decoherence on the performance of the 
optimal Toffoli gate has been explored (Figure 3), and 
we interpret the corresponding result as follows: when 
T 3> 0, the effect of thermal noise becomes less dom¬ 
inant, and the fidelity of a quantum gate is almost en¬ 
tirely contributed by the intrinsic fidelity. Here we have 
been able to design a fast and optimal pulse for a Tof¬ 
foli gate using our quantum control approach for which 
the intrinsic fidelity is so high (~ 99.99%), that the fi¬ 
delity (F ~ 99.9%) is significantly contributed by the 
decoherence-induced noise, which is also very small (com¬ 
pared to previous realizations) for the state-of-the-art su¬ 
perconducting atoms with T ~ 20 — 60 /rs [3] . 

In summary, we have devised a powerful quantum con¬ 
trol scheme, named SuSSADE, to design a fast and high- 
fidelity single-shot Toffoli gate for a scalable chain of 
nearest-neighbor-coupled three-transmon system. The 
time required for the Toffoli operation is comparable 
with the timescale of two-qubit avoided-crossing-based 
CZ gate, which is the key advantage our quantum-control 
approach proffers compared to decomposition-based ap¬ 
proaches requiring many such two-qubit gates to imple- 
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merit a single CCZ operation. Our three-transmon sys¬ 
tem serves as a module for all ID and 2D quantum com¬ 
puting architectures [34], and therefore one can realize 
our scheme in a large-scale multi-qubit architecture, if 
the undesired couplings are turned off [35]. Our ap¬ 
proach demonstrates the efficacy of SuSSADE for design¬ 
ing quantum gates as well as yielding the concrete exam¬ 
ple of a three-qubit gate required for scalable quantum- 
error-correction. 
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SUPPLEMENTAL MATERIAL FOR 
“HIGH-FIDELITY SINGLE-SHOT TOFFOLI 
GATE VIA QUANTUM CONTROL” 

Subspace-Selective Self-Adaptive Differential 
Evolution (SnSSADE) 

Differential Evolution (DE) was initially introduced by 
Storn et al. [44] as a global optimization approach and is 
currently known to be the most efficient of all the evolu¬ 
tionary algorithms [45]. In its basic form, DE performs 
three operations namely: mutation, crossover and selec¬ 
tion. These three operations cooperatively evolve the 
IV-dimensional candidate solutions Ci, f C {1,..., P} to¬ 
ward their optimal position with N and P being the num¬ 
ber of optimization parameters in the search space and 
population size respectively. 

The elements of Ci are piecewise-constant functions, 
where each element is fixed over the time interval Ar = 
Q/N, with 0 being the total time needed to perform the 


unitary operation. DE mutates the candidate solution as 
follows: 


Di — Ci.^ -\- fi {Ci^ — ) (4) 

where ii, * 2 , * 3 ) * C [1;^] ^re integers and mutually dis¬ 
tinct. The quantity jr G [0, 2] is a uniformly distributed 
random number, which defines the step size by which 
DE explores the search space, and Di is the donor vector 
resulting from the mutation operation (4). 

The next step is to switch the elements Di{j) of donor 
vector and each candidate solution Ci via the crossover 
operation 


W) 


Di{j), if X < C 
Ci{j), otherwise 


( 5 ) 


with X G [0,1] is a random number chosen from a uniform 
distribution and ^ defines the crossover rate. The final 
step is the selection 


5 ' = 



if F{C,) < f{X,) 
otherwise 


( 6 ) 


with S'i being the offspring of Ci for the next genera¬ 
tion and f{Xi) being the objective function, which is the 
intrinsic fidelity P for our case. 

Here the standard version of DE was not able to deliver 
a fidelity better than 95% for the Toffoli gate. Thus, 
we enhance the standard DE algorithm in order to find 
the optimal solution for problems with high dimension. 
The first step towards this enhancement is to find the 
optimal value of algorithmic parameters, namely /r and 
^ [40] . One approach is to try DE with many trial guesses 
for the algorithmic parameters to find the best set for /j, 
and leading to the best optimal solution for a specihc 
problem under consideration. However, this trial-and- 
error approach becomes intractable (as it is in our case) 
when the computational cost of evaluating the objective 
function is expensive, and running DE for a long time 
requires an excessive amount of computational resources. 

As an alternative approach, we use the self-adaptive 
DE [36] , which self-adaptively updates the value of /i and 
^ at each generation C using 


and 




lH+ri.iiu ifr2<«:i 
\ii^G otherwise 



if r4 < 
otherwise. 


( 7 ) 


( 8 ) 


thereby producing the new algorithmic parameters in a 
new parent vector (candidate solution in each genera¬ 
tion). The quantities r^, with j G 1,2, 3,4, are uniform 
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random numbers generated from [0,1], and ki and K 2 
are the probabilities to adjust the algorithmic parame¬ 
ters. The quantities ki, K 2 , fJ-u and jii are assigned to 
fixed values: 0.1, 0.1, 0.1 and 0.9 respectively. Whereas 
using self-adaptive DE by itself improves the obtained 
fidelity, the resulted fidelity still remains under the req¬ 
uisite threshold 99.9%) due to the high-dimensionality 
problem. 

The next step to enhance the standard DE is to 
combine self-adaptive DE with cooperative coevolution 
(DECC-II) [33], designed for high-dimensional optimiza¬ 
tion problems. Even the combination of self-adaptive DE 
with DECC-II in their original forms failed to yield a bet¬ 
ter fidelity so modified this combination. Here we first 
describe the DECC-II followed by our enhanced version 
of DE. 

DECC-II works by decomposing the iti-dimensional 
candidate solutions to m-dimensional subspaces and op¬ 
timizes each subspace for s cycles at each generation, 
whereas other subspaces remain unchanged. DECC-II 
uses Self-Adaptive Differential Evolution with Neighbor¬ 
hood Search (NSDE) [46] to optimize each individual 
subspace. The number s of cycles and the dimension m 
of the subspace should be assigned by the user, and se¬ 
lecting the optimal value is a computationally expensive 
task. We found the results to be largely influenced by 
the choice of s and m, with none of them satisfying our 
requirement of having a threshold-fidelity Toffoli gate. 

Inspired by the DECC-II algorithms, we set s = 1 and 
choose the dimension of subspace m randomly from {1, 
2, 3, 4, 5}. Now, we no longer need to look for the op¬ 
timal values of s and m for the optimization. However 
this new approach has an extremely slow convergence, as 
in each generation only a small part of the candidate so¬ 
lutions are being selected for the optimization. In order 
to speed up the convergence, our algorithm randomly 
switches breeding between the subspace and the whole 
space according to the value of an input switch parame¬ 
ter S S [0,1], such that a uniformly distributed random 
number rj € [0,1] at generation j restricts breeding to 
the subspace, if rj < S, and breeds in the whole space 
otherwise. As our algorithm selects an m-dimensional 
subspace at each generation and self-adaptively evolves 
the algorithmic parameters, we call it Subspace-Selective 
Self-Adaptive DE (SuSSADE). For our purpose, we ob¬ 
serve that choosing m = 1 suffices, which signihes that 
the selected subspace is trivial. We refer to this one¬ 
dimensional extreme case as IDSUSSADE. 


Phase compensation 

In order to design the optimal pulse for the CCZ gate, 
we employ another strategy, called “Phase Compensa¬ 
tion”, which is frequently used in designing quantum 
gates with superconducting atoms [16, 21]. Any arbitrary 


rotation about the z-axis is trivial for a superconducting 
qubit and is usually performed via qubit frequency ex¬ 
cursions. Therefore, we assume that any phase acquired 
by |100), |010), and |001) states can be nullified by post 
z-rotations. 

Whereas an ideal CCZ gate is given by (in tensor 
product basis as dehned in the main text), CCZideai = 
diag(l, 1,1,1,1,1,1, —1). Under the freedom that 
some additional phases can be compensated with post¬ 
rotations about z-axis, it is sufficient for our purpose if 
we dehne our target gate as 


CCZ := diag (l, 

( 9 ) 


_^i{6i+0'2+63)\ 


with 9j being any arbitrary phase acquired by the j 
transmon under the optimal pulse. 


th 


Decoherence 

In order to estimate the effect of decoherence, we model 
each superconducting transmon as a damped harmonic 
oscillator and then compute the time evolution of the 
density matrix under amplitude and phase damping [2, 
22]. We employ Kraus’s operator-sum representation to 
represent the completely-positive trace-preserving maps 
for these damping processes. If p{t) is the density matrix 
of the system at any given time t, then, under a specihc 
dissipative process, 


p{t) = 2/cWp(0)<£lW> (10) 

fc =0 

with £fe’s the Kraus matrices for the given process satis- 
fying = 1 at each time instant. 

We now discuss construction of Kraus matrices for am¬ 
plitude and phase damping of a single transmon. Once 
we know the Kraus matrices for each transmon, the set 
of Kraus matrices for the entire three-transmon system 
can then be constructed as all possible tensor products 
of those single-transmon Kraus matrices, assuming that 
decoherence affects each superconducting atom indepen¬ 
dently. 

Amplitude damping: The Kraus matrices for ampli¬ 
tude damping of a single superconducting atom are given 
by (assuming 4 energy levels for each atom) [22] 

j—l V 

( 11 ) 

with I e {0,1,2,3} denoting the indices for the 
Kraus matrices. It is straightforward to verify that 
ul(t)i?/(t) = 1 at all time. The exponential 
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amplitude-damping factor exp(—t/Ti) has relaxation 
time Ti- 

Phase damping: The Kraus matrices for phase damp¬ 
ing of a single superconducting atom are given by (as¬ 
suming 4 energy levels for each atom) [22], 


3 


Ez(t) = y^exp 





n 


\j) 01 : 


( 12 ) 


with I again denoting Kraus-matrix indices. Here T 2 
denotes the dephasing time. For phase damping, 
= 1 at all times is straightforward to 
check. Such a completeness relation cannot be satisfied 
if we truncate the infinite series, thereby constraining the 
Kraus-matrix index I to a finite upper bound. We must 
impose an upper limit for I in order to ensure a finite 
number of Kraus matrices. Such a truncation is valid as 
long as our total gate time is much shorter than the de¬ 
phasing time T 2 (implying t/T 2 ^ 1), which is true for 
our case. In our calculation here, we assume 0 < ^ < 3 for 
phase damping as well, neglecting all higher order terms 
for t/T 2 . 


Proposed experimental realization 

In this section we discuss the calibration procedure to 
suppress the effect of residual noise in an experiment, 
where a fast and high-fidelity Toffoli gate is demon¬ 
strated. In a quantum control procedure, the effect of 
noise is to distort the external pulses e.g. time-dependent 
qubit frequencies, which are designed optimally through 
a control procedure. In an ideal experiment with no noise 
acting on the quantum system. Fig. 1(a) shows our op¬ 
timal pulses in the absence of distortion. However this 
idealized assumption is not valid in a real experimental 
procedure when noise distort the designed pulses. There¬ 
fore a calibration procedure must take the distortions into 
account for optimally designed square pulses. 

Up to first order, distortion occurs in an experiment 
whenever a square pulse is generated and passed through 
the superconducting control electronics. Therefore, by 
the time that the square pulse reaches the artificial atom, 
the pulse becomes smoothed. Mathematically, we model 


first-order distortion by convolving the square pulse with 
a Gaussian function, which yields an error-function (erf) 
pulse [16, 21]. Leading-order distortion has already been 
taken into account in designing our optimal pulses as 
shown in Fig. 1(b) in the original manuscript. 
Calibration is required for higher orders and can be 
achieved through the well-known Closed-Loop Learning 
Control (CILC) technique [37, 42]. CILC has been suc¬ 
cessfully applied to tasks such as discriminating simi¬ 
lar molecules [39], ionization [43], molecular isomeriza¬ 
tion [38], and coherent quantum control of two-photon 
transitions [41]. 


Noise 

I 



FIG. 4: (color online) A schematic view of the Closed-Loop 
Learning Control (CILC) technique which uses SuSSADE as 
a quantum control scheme. Based on this control procedure, 
the optimal pulses which are generated using SuSSADE fed 
into a noisy quantum system, if the output met the target 
fidelity, the procedure aborts otherwise a new set of control 
pulses is chosen using SuSSADE, and the iteration continues 
until the target is met. 


CILC is an iterative technique for searching the optimal 
solution in a quantum control landscape. The idea is 
to start with a quantum control scheme (SuSSADE for 
our case) to solve an ideal, physically realistic model of 
the quantum system and then test the obtained optimal 
solution experimentally. If the experimentally measured 
objective function (fidelity in our case) does not satisfy 
the target, the control pulse should be calibrated again 
with the same quantum control algorithm (SuSSADE for 
our case) in order to obtain a new set of optimal solutions, 
and this process continues until the target is met. A 
schematic view of CILC is depicted in Eig. 4. 









